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Several f[R) modified gravity models have been proposed which realize the correct cosmological 
evolution and satisfy solar system and laboratory tests. Although nonrelativistic stellar configura- 
tions can be constructed, we argue that relativistic stars cannot be present in such f{R) theories. 
This problem appears due to the dynamics of the effective scalar degree of freedom in the strong 
gravity regime. Our claim thus raises doubts on the viability of f{R) models. 
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I. INTRODUCTION 

The current accelerated expansion of the Universe is 
one of the deepest mysteries in cosmology T] . This accel- 
eration may be due to some unknown energy-momentum 
component having the equation of state p/ p w —1, or 
may be due to a modification of general relativity. In 
this paper, we are interested in the latter possibility. The 
simplest phenomenological way of modifying gravity is to 
consider a gravitational action described by a function of 
the Ricci scalar, f{R), instead of the Einstein-Hilbert 
action. An early attempt is found, e.g., in [3|, where a 
modification like f{R) — R-{- R^/fJ.^ was used to explain 
the accelerated expansion in the early Universe. More 
recently, f{R) modified gravity theories are often consid- 
ered as a possible origin of the current acceleration of the 
Universe [3|. 

Any modified theories of gravity must account for the 
late-time cosmology which is well established by observa- 
tions, and at the same time must be consistent with solar 
system and laboratory tests of gravity. However, since 
f{R) gravity has the equivalent description in terms of 
the Brans-Dicke theory with the Brans-Dicke parameter 
a; = [4| , naively constructed models would result in vi- 
olation of the above requirements [a, S Q- ^^''^ example, 
the original proposal of \^ employs f{R) = R — p,^/R, 
which admits an acceleratedly expanding solution even in 
the absence of a dark energy component. In order to ac- 
commodate a late time acceleration, however, one must 
introduce a very low energy scale, /i ~ Hq (the present 
Hubble scale), leading to a very light scalar field, which 
predicts the parametrized post-Newtonian (PPN) param- 
eter 7 = (1 -I- a;)/(2 -\- uj) = 1/2. Obviously, this result 
contradicts the observational constraint I7— 1| < 10^** [9] 

To circumvent this difficulty, it is important to no- 
tice that the presence of matter may affect the dy- 
namics of the extra scalar degree of freedom. The key 
idea is essentially the same as that of the "chameleon" 
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model [10|, [ill, II4I, in which the effective mass of the 
scalar field depends on the local matter density. In 
particular, the scalar field is very light for the cos- 
mological density and is heavy for the solar system 
density, though the actual mechanism is slightly more 
complicated. The most successful class of f{R) mod- 
els O, H E [H, [13, [ll E H [in incorporates this 
chameleon mechanism to evade local gravity tests. The 
experimental and observational consequences of this kind 
of f(R) models are found in Refs. [23, IH ll^. (See 
also i25|.) 



In this paper, we consider the strong gravity regime 
of the carefully constructed models of [T^, [20, [2l| ■ The 
strong gravity aspects of f{R) theories have not been ex- 
plored so much before. Recently, Frolov suggested that 
such f{K) models generically suffer from the problem of 
curvature singularities which can be easily accessed by 
the field dynamics in the presence of matter [26|. In 
other words, a curvature singularity may be caused not 
by diverging gravitational potential depth, |$| = cx), but 
rather by a slightly large gravitational field, |$| < 1/2. 
This motivates us to study relativistic stars in the con- 
text of f{R) gravity. Spherically symmetric stars in f(R) 
gravity have been investigated so far in [27|, [28|, [29|, [30| . 
(We confine ourselves to a metric theory of f{R) gravity. 
Using the Palatini formalism, polytropic stars have been 
studied in |3l|.) We shall show, both analytically and 
numerically, that stellar solutions with relatively strong 
gravitational fields cannot be constructed. Using the spe- 
cific example of relativistic stars, we clarify how the sin- 
gularity problem arises in the strong gravity regime of 
f{R) theories. The singularity problem was also identi- 
fied in [32] in a cosmological setting. 

This paper is organized as follows. In the next section 
we derive equations of motion for f{R) modified grav- 
ity, and define the specific model we study. In Sec. IIIIl 
we reinterpret the problem of finding the desired stellar 
configuration as the problem of the particle motion in 
classical mechanics. We give some analytic arguments 
in Sec. IIVI and then we present our numerical results in 
Sec. |Vl Finally, we draw our conclusions in Sec. IVII 



II. PRELIMINARIES 



A. f{R) gravity 



The action we consider is given by 

fiR) , 



S = d xy/—g 



IGttG 



C„ 



(1) 



where f{R) is a function of the Ricci scalar R, and £m is 
the Lagrangian of matter fields. Variation with respect 
to metric leads to the field equations 

fRR^.. - ^t^/ufR + f n/fl, - ]:f] g^. = 8^GT,,,, (2) 



where Jr :— df/dR and T^i/ := 
The trace of Eq. ([2]) reduces to 



-2SC^/Sg^''' + g^,C, 



OfR 



SttG 



1 



T+-i2f-fRR). 



(3) 



We now introduce an effective scalar degree of freedom 
by defining x := f^. Inverting this relation, the Ricci 
scalar can be expressed in terms of x- R — Qix)- Thus, 
Eqs. ([2]) and ([3]) are equivalently rewritten as 



Ox 



SttG, 



2x^ dV 



-T+-^^ 



where the effective potential V is given by 



and 



v{x) 



dV 



^ ixQix) - fiQix.))] , 



rf^-^[2/(Q(x))-xQ(x)] 



(5) 



(6) 



(7) 



Eqs. (j4|) and ([5|) are equivalent to the Jordan frame 
equations of motion in the Brans-Dicke theory with uj = 
0, if we ignore the potential term V{x)- One can move 
to the Einstein frame by performing the conformal trans- 
formation g^^ = XQtif with x = exp(A/l67rG/3 0), where 
(f) is the canonical scalar field. The potential for (f) is then 
given by V{x{4'))- Although the Einstein frame equa- 
tions of motion are sometimes convenient, we shall not 
use this and work in the Jordan frame directly through- 
out the paper, so as to avoid confusion concerning the 
coupling between matter fields and the effective scalar 
degree of freedom. 

B. The model 

In order to be explicit, we take f{R) in the following 
form [11: 



f{R)^R + XRo 



R^ 
Ri 



(8) 






1.125 

1.1 

1.075 

1.05 

1.025 

1 

0.975 




0.7 0.75 0.8 0.85 0.9 0.95 1 

X 



FIG. 1: The effective potential V{x) for Starobinsky's f{R) 
model with n — 1 and xi — 3.6 (A ~ 2). x is the effective 
scalar degree of freedom defined by x •= df/dR. 



where n, A(> 0), and Ro{> 0) are parameters. This 
model is carefully constructed so that it gives viable cos- 
mology and satisfies solar system and laboratory tests. 
Suppose that the de Sitter solution in this theory is ex- 
pressed as R = Ri = constant — xiRq. We may define 
the "cosmological constant" as A^s := i?i/4. Since the 
de Sitter solution follows from dV/dx = 0, one has 



A^i- 



xi{i + xi) 



2\n+l 



(9) 



2(l + a;^)"+i- 1- {n+l)xl' 

Thus, we may use xi as a model parameter instead of 
A. We will take as xi the maximal root of Eq. ^ for a 
given A, because it corresponds to the de Sitter minimum 
of the potential. 

The scalar field x is written in terms of R as 



X 



l-2nX 



R 

Rn 
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R^ 
Rl 



(10) 



One sees that x ^ 1 as i? — » ±oo and R^ Q. Note that 
curvature singularities, R = ±cx), is mapped to the finite 
value of X = 1 • The value of x at the de Sitter minimum 
is given by 



XI = 1 



(1-f a;f)"+i - 1- (n + l)a 



(11) 



A typical form of the potential V{x) around the de Sit- 
ter minimum is shown in Fig. [TJ In fact, the potential 
is a multivalued function of X: and its shape is compli- 
cated, probably even pathological, away from the plotted 
region. However, since our discussion here focuses on the 
behavior of x around the de Sitter minimum, there is no 
difficulty with such a complicated potential. 

Though we shall confine ourselves to the specific model 
defined by Eq. ([5]), the result will apply to other f{R) 
models as well. In particular, the models proposed by 
Hu and Sawicki \2(^ and^by Appleby and Battye [2l| fall 
into the same class as |19| in the sense that the potential 
around the de Sitter minimum has the same structure. 



III. SPHERICALLY SYMMETRIC STARS IN 
f{R) GRAVITY 

A. Basic equations 

To study stellar configurations in /(i?) gravity, we take 
the ansatz of a spherically symmetric and static metric; 



ds^ = -N{r)dt^ + 



1 



B{t) 



dr^ + r^ ide'^ + sin^ OdLp^) .(12) 



The energy-momentum tensor of matter fields is given by 

t; ^di&g{- p,p,p,p). (13) 

From the energy-momentum conservation, V^TJ^ = 0, 
we obtain 



N' 



(14) 



Here and hereafter a prime denotes differentiation with 
respect to r. The (tt) and (rr) components of the field 
equations (|4]) yield, respectively, 



4 (-1 + S + rB') = -SnGp - x^V 



-B 



2 B'\ , 



(15) 



4 f -1 + B + rB^ ] = SnGp - x^V 



N 



-«I; + S^U' («) 



Note that using the scaling freedom of the t coordinate, 
we set 7V(0) = 1. From Eqs. ^-^, we obtain 



3^2 

B2 -h 2N2 



-8nGpc - XcVc - 3C2, (19) 

inGpc - XcVc - 2C2, (20) 



8ttG, 



2x1 



3C2 = -^i-Pc + 3pc) + -^V^^, (21) 

where G := G/xc is the effective gravitational constant, 
Vc :— V{xc), and V^^ = dV/dx\x=xc- These three equa- 
tions are rearranged to give 



B2 

N2 



^■^G Xc 

-g- (2pC + 3PC) - -J 



2x| 
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Vc-=^Vx^, (22) 



—^{2PC + 3PC, g 



^V - —V 
" 9 ^=' 



SvrG 



2X^ 



^2 = -^i-Pc + Spc) + -^Vx^. 



(23) 
(24) 



Then, p2 is derived from the conservation equation: 

P2+N2{pc+Pc)^0- (25) 

The Ricci scalar is given hy R = Re + 0{r^) with Re = 
-6(^2 + N2) near r = 0. 

If the energy density is constant inside the star, p = po, 
Eq. (fT4)) immediately gives 



A^(r) = 



PO+Pc 



Po+p(»-) 



(26) 



In the rest of the paper, we focus on constant density 
stars for simplicity. 



The equation of motion for x [Eq. ^] gives 
B 



X 



2 N' B' 

r^2N^2B 
SnG 



X 



2x'dV 
3 (-^ + '^) + — d^- ^^'^ 

We will not integrate the angular components of the 
field equations. Instead, we will use them to check the 
accuracy of our numerical results, because those are de- 
rived from other equations via the Bianchi identity. 

To specify the boundary conditions at the center of a 
star, assuming the regularity, we expand the variables in 
the power series of r as 



N{r) = 1 + N2r'^ + ..., B{r) ^ 1 + Bar^ -|- 



(18) 



p(r) 



, P2 2 
Pc + ^r 



•, P{r) 



, P2 2 
Pc + ^r 



where Xci Pc and Pc are the central values of the scalar 
field, the energy density and the pressure, respectively. 



B. A classical mechanics picture 

Given poi Pc, and Xc (or, equivalently. Re), Eqs. p^ 
([TT]) can be integrated outwards from the center to the 
surface of the star, r = TZ, which is defined by p{TZ) — 0. 
Then, one integrates the vacuum field equations from 
the star surface to sufhciently large r, finding the exte- 
rior profile of the metric and the scalar field. Unlike in 
general relativity, we have the extra scalar degree of free- 
dom corresponding to the choice of Xc- However, not all 
values of Xc can lead to a physically reasonable solution 
inside and outside a star. The desired solution is such 
that x(r) — > xi as r — > CX3, i.e., the asymptotically de 
Sitter solution with the cosmological constant Acff. 

We can formulate the problem of finding the physical 
configuration of stars as follows. Eq. p?|) can be written 
as 



d^ 
dr"^ 



2 d ^ 
r dr' 



dU 
dx 



^, 



where 



dU 
dx 



UrQ - 2/] 



2x^dV 
3 dx 



(27) 
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FIG. 2: The (inverted) potential U(x) for Starobinsky's f{R) 
model with A = 2 and n — 1. The point A corresponds to a 
curvature singularity {R — -|-cxd), and the point B is the de 
Sitter extremum. (See also Fig. 1 of Ref. :26i|.) 
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FIG. 3: Motion of a particle near the de Sitter extremum 
of f/(x)- The particle feels the force T (< 0) which arises 
from the trace of the energy-momentum tensor of the matter, 



and 



^:=-^(p-3p). 



(29) 



Here we neglect the effect of the metric for the moment 
to comprehend the essential point. [Later we will solve 
the full set of equations P^ - (|17p numerically.] 

Regarding r as a time coordinate, we find that Eq. ([27]) 
is a "dynamical" equation describing the motion of a par- 
ticle in the potential IJ under the time-dependent force T . 
The second term in the left hand side of Eq. (|27p repre- 
sents frictional force, which may also affect the dynamics 
in some cases. The potential for the particle t/(x), de- 
fined by Eq. ((28|) , is different from the inverted potential 
— y(x). However, the structure of U around x = Xi is 
quite similar to —V, as is shown in Figs. [2] and [S] the 
point A (x = 1) corresponds to a curvature singularity, 
i? — > CX3, and the point B (x = Xi) is the de Sitter ex- 
tremum. 



In Ref. [201, Frolov similarly introduced the potential 
— C/(x) and the force term which is essentially given by 
the trace of the energy- momentum tensor. For the pur- 
pose of solving for the radial profile, it is more convenient 
to consider the inverted potential — (— f/) = ?7, as in 
the cases of bubble nucleation [33| and of the chameleon 
model [10] 

Suppose that the initial position of the particle, Xci lies 
between points A and B. The particle starts at rest since 
;;;^'|^^Q = 0. The force term T depends on the matter 
configuration inside the star and plays a crucial role in 
this problem. When p > 3p, we have J^ < 0. Since 
the pressure becomes smaller for larger r, the force [JFI 
is stronger near the surface than in the central region of 
the star. The force vanishes for r > 7?. (if one assumes 
the vacuum exterior). 

For fixed po and pc, the behavior of the particle de- 
pends on its initial position. If Xc is sufficiently close to 
1, the slope of the potential is bigger than the force term 
J^ initially, so that x "^iH rapidly roll down to the cur- 
vature singularity, x = 1- Let Xs be the minimum value 
of Xc for which this occurs. Specifically, Xs is found by 
solving the equation [T — dU/dxlxs =0. If Xc < Xs, 
the particle initially climbs up the potential hill under 
the force ^ < 0. The force vanishes a.t r = TZ, but the 
particle keeps climbing up the potential for the moment. 
For r > TZ, the scenario one can easily deduce is as fol- 
lows: The particle cannot reach the top of the potential 
hill and turns around at some r, ending up with the cur- 
vature singularity x = 1, or goes through the top and 
rolls down to the left, depending on the initial position 
Xc [35l |. In this case, one finds a critical value Xcrit be- 
tween the "turn-around" and "rolling-down" solutions. 
By fine-tuning the initial position so that Xc = Xcrit, 
we can realize the asymptotically de Sitter solution for 
which X ~* Xi as r -^ cx). Thus, the problem reduces to 
a boundary value problem. 

However, there is another possibility that the particle 
inevitably overshoots the potential even for the possi- 
ble maximum value of Xc (i-e., Xs)- In this case, one 
cannot obtain the desired solution: the particle rolls 
down into the curvature singularity right after it starts 
to move, or overshoots the potential. Since the gravi- 
tational potential produced by a star is proportional to 
{p()R^)/TZ = poTZ^, a stronger gravitational field implies 
stronger force and/or a longer period during which the 
force term survives effectively. As a result, the fine-tuned 
initial location Xcrit goes toward the right as the star ac- 
commodates a larger gravitational potential, and Xcrit 
will eventually reach the point Xs- Therefore, it is ex- 
pected that there is a maximum value of the gravitational 
potential for a star to exist. 



IV. ANALYTIC ARGUMENT 

Before solving the full set of equations (|14 p - p7|) nu- 
merically, in this section we shall provide some analytic 



arguments. The analysis with approximate solutions will 
help to understand our numerical results presented in 
Sec. El 

First let us consider the interior of a star: r < TZ. To 
give a tractable argument, we assume that 

|S-1|, |iV-l|<l, \B'/B\, \N'/N\-^r-\ (30) 



Then, the solution to Eq. (fT7| which is regular at the 
center is given by 
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(31) 



(32) 



Here we have ignored the pressure p relative to poj a-nd 
made a rough approximation x^^x ~ X^^xc- Using 
Eqs. (US]) and ^, we obtain 



^^^_2G(M_mZ3) 

7^3 



(33) 



where M := AnpoTZ^ /3. Here we have neglected the "cos- 
mological constant" x^^- (We remind the reader that 

G:=G/xc.) 

To derive the exterior solution, we approximate 



3 dx 



x-xi 



(34) 



^ .— {2xi/'i)d^V/dx\xi^ and analyze the be- 



where A^^ : 

havior of x around xi- The exterior solution to Eq. (|17|) 

is found to be 



X-Xi+C' 



-(r-n)/\^ 



(35) 



Matching this to the interior solution (pi]) at r = TZ, we 
obtain 



C 



2G/i 



Xi 



G^ 

-1/2n 



(36) 

(37) 



where we have used A,^ ~ ©(A^'^) :$> TZ. Given po, 
Eq. (P7|) determines Xc- Solving Eq. P^ and matching 
the solution to the interior one (j33p at the surface of the 
star, we find 



B~ 1- 



2G{M - n/3) 
r 

Then, Eq. (fT6l) implies that 

2G'(Af + ^/3) 



N-Nr. 



1 



(38) 



(39) 



where the asymptotic value Noo is not unity because of 
our boundary condition at the center, but it can be set 
unity by rescaling of the time coordinate. From Eqs. (|38|) 
and (|39p . the PPN parameter 7 is found to be 



7 



MI - n 
3M 



fi 



(40) 



Now let us define the gravitational potential evaluated 
at the surface of the star: 



$ 



G{M + ^l/3) 
TZ 



In terms of this, Eq. (|37p can be written as 



B/i 



* 3M + /i 



with A := 



Xc-Xi 

Xc 



(41) 



(42) 



The "thin-shell" condition is given by A ^ $ 10]. This 
is equivalent to /i <C M, which leads to C ^ GM and 
indeed suppresses the deviation of x(^) from xi outside 
the star. This situation is realized if AnGpo 



xiv^. 



0. 



On the other hand, if inGpo ^ xt^xa ; and hence p ~ M, 
the thin-shell condition does not hold. It is easy to see 
that 7 ~ 1/2. In this "thick-shell" case, Xc must be 
not too far from xi so as not to fall into the curvature 
singularity, x — ^- Therefore, we require that Xc < 1- 
This condition together with Eq. ([57)1 gives the bound 



$ < $„ 



Ki-xi), 



(43) 



i.e., stars with <I> > $,nax cannot exist. As shown later by 
numerical solutions, the thin-shell condition is violated as 
long as the exterior is vacuum. 

Specifically, for n = 1 one needs xi > \/3 in order 
to have a de Sitter minimum. This leads to 1 — xi < 
1/3 ($max < 4/9). For xi = 3.6 (A = 2.088), we have 
1 - xi = 0.0771 6 ($max = 0.1029). For n = 2, one needs 
xi > \/yi3-2, giving 1-Xi < 0.2705 ($max < 0.3606). 



V. NUMERICAL RESULTS 



A. Stars in f{R) gravity 

We numerically integrate Eqs. p^ ^ P?)) outwards from 
the center, imposing the appropriate boundary condi- 
tions given by Eq. HH]) with Eqs. dHjl-dinD. We recon- 
struct R — R{r) from the metric and compare it with the 
solution for the scalar degree of freedom, R = Q{x{^))i 
to make sure that numerical errors are sufficiently small. 
The angular component of the field equations is also used 
for the same purpose. 

The sets of the model parameters we use are: (i) rt = 1 
and xi = 3.6 (A = 2.088), for which xi = 0.9228; (ii) 
n = 2 and a;i = 3.6 (A = 1.827), for which xi = 0.9903. 
For numerical solutions we shall take AirGpo = lO^Aoff. 
This is not a realistic value, e.g., for neutron stars, but 
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FIG. 4: Metric for a nonrelativistic star. Parameters are given 
by n = 1, xi = 3.6, A-uGpo = lO'^Acff , and pc = 5 x lQ~^po. 
The central value of the Ricci scalar is tuned to be R^ — 
3.462 X 10~^ X SirGpo- The radial coordinate is normalized 
by the radius of the star 7?.. 
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FIG. 5: Numerical solutions of the Ricci scalar and x &* 
a nonrelativistic star. Parameters are the same as those in 
Fig.|4] Dashed line is a plot of the analytic approximation (I3ip 
and (|35|) . Since R ^ xiRq and x ^ Xi ^^ ^ ~^ °Oj Ihis is the 
desired solution with asymptotically de Sitter geometry. 



we do not have to be concerned about this point because 
the properties of stellar solutions which we are interested 
in are basically characterized by the gravitational poten- 
tial rather than the energy density itself, provided that 
Gpo 3> Aoff. This is expected from the analytic result, 
and we have confirmed that it is indeed true for our nu- 
merical solutions. 



1. n = 1 and xi = 3.6 

Taking 47rGpo = lO^Aofi and Pc = 10~^/9o, we find 
an asymptotically de Sitter solution of a star for Re — 
2.001 X 10~^ X SnGpo- The solution agrees well with the 
analytic approximation in Sec. IIVI 

A numerical calculation has been performed also in 



FIG. 6: The pressure profile p{r) of a nonrelativistic star. 
Parameters are the same as those in Fig. |4l 



the case of AnGpo — lO^Acff and pc = 5 x 10 ^po- We 
find a solution by tuning Re = 3.462 x 10^^ x SnGpo 
{Xc = 0.9836), for which 



47r 



Gpo7^' = 



GM 



0.06687. 



(44) 



Our numerical result is shown in Figs.[4]-[6l One can see 
from Fig.Othat R — > xiRq and x — ^ xi as r — > oo. Since 
A — 0.06176 is almost the same as the gravitational po- 
tential (UH), the thin-shell does not form. A numerical 
fitting leads to the approximate expression for the exte- 
rior metric: 



N 
B 



1- 



l-2ci^-|Aeffr2 
r 3 



2c3- - ^Aeffr^, 
r 3 



(45) 
(46) 



with 



N^ = 1.332, ci = 0.08716, ca = 0.9973, 

C3 = 0.04747, Ci = 0.9993. (47) 

The PPN parameter turns out to be 7 ~ cs/ci ~ 0.5446, 
which is close to 1/2, as expected. 

To explore stars with larger GM/TZ, we have tried to 
find numerical solutions for AirGpo = lO^Aoft and Pc = 
0.1 X pq. a solution which is regular inside the star is 
obtained, e.g., for Re = 0.7000 x SirGpo, but this is the 
rolling-down solution (see Fig. [7]) . Hence it is unphysical 
{N{r) -^ and B{r) ^ 00 as r —^ 00), though the 
gravitational potential is as "large" as 



GM 



0.1203. 



(48) 



Taking a slightly larger value of the central curvature, 
Re — 0.7001 X SnGpo, the Ricci scalar rapidly diverges 
inside the star. This corresponds to the case in which x 
starts to move to the right toward the curvature singular- 
ity, X = 1- Note that since x^V « R/2, Re = O.TxSnGpo 
implies 47rG(po - 5Pc) - xl^xa ~ 0. 
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FIG. 7; A rolling-down solution for a would-be relativistic 
star. Parameters are given by n = 1, ii = 3.6, 4nGpo — 
lO^Acff, and pc = 0.1 x po. The central value of the Ricci 
scalar is Re — 0.7000x87rGpo. The solution clearly overshoots 
the de Sitter extremum, X = Xi- 



2. n — 2 and xi 



3.6 



For 4:TtGpo = lO^Aeff and pc = 5 x 10 '^po, we find the 
stellar configuration with de Sitter asymptotic behavior 
by tuning R^ = 2.035 x IQ-*^ x SnGpo (Xc = 0.9911). In 
tfiis case, GM/TZ = 7.491 x 10"". The thin-shell condi- 
tion does not hold (A = 7.448 x 10"''), and 7 = 0.5005. 

Looking for stars with stronger gravitational fields, we 
take AnGpo = lO'^Acff and pc = I0~'^po. When R^ == 
0.9700 X SirGpo, a rolling-down solution is found with 
GM/n = 0.01465, while a slightly larger value R^ ^ 
0.9701 X SirGpo leads to a curvature singularity R ^ 00 
inside the star. 



From the above numerical analysis, we conclude that 
stars with strong gravitational fields cannot be present 
in this class of f{R) theories. 



B. Stars in surrounding media 

So far we have examined the case of the vacuum exte- 
rior. Although we can indeed construct a stellar configu- 
ration provided that gravity is weak, the exterior metric 
is not given by the de Sitter-Schwarzschild solution in 
general relativity and the PPN parameter is found to be 
7 ~ 1/2. This simply reflects the fact that the thin-shell 
condition is violated in the present case. However, we 
can make the chameleon mechanism effective by taking 
into account the effect of surrounding media (e.g., dark 
matter). Let us make a brief comment on this point. 

In a "realistic" situation, exterior matter is present 
around a star, giving rise to the force term J-Qxt — 
— 87rG/9cxt/3 there. Then, to obtain a viable stellar con- 
figuration, one has to take a shot at the point x* satisfy- 



ing the equation jT^xt — dU/dx = rather than the top of 
the potential hill, xi- Since x* > Xij the difference in x 
between inside and outside the star becomes smaller, and 
hence it is easier to satisfy the thin-shell condition. In- 
deed, the chameleon mechanism has been shown to work 
in the f{R) model which is very similar to the current 
one, reproducing 7 ~ 1 in the solar vicinity [20|. 

The above argument, however, only applies to stars 
with weak gravitational fields. When gravity is strong, 
for any initial condition of Xc which avoids the curvature 
singularity, the scalar field x inevitably overshoots the 
potential. It is clear from this fact that one cannot stop 
X at any value of x* > Xi- Therefore, a relativistic star 
cannot be present even with surrounding medium. 



VI. CONCLUSIONS 

In this paper, we have studied the strong gravity aspect 
of /(i?) modified gravity models that reproduce the con- 
ventional cosmological evolution and evade solar system 
and laboratory tests [1^, [20, [21| . It is known that f{R) 
theories can be recasted simply in the Brans-Dicke the- 
ory with uj — 0, but the potential for the effective scalar 
degree of freedom may play a complicated and nontrivial 
role. Moreover, the presence of matter may affect dynam- 
ics of the scalar field, possibly mimicking the chameleon 
model [10]. 

We have explored uniform density, spherically symmet- 
ric stars and their exterior geometry in the f{R) model 
of [19]. The main result of the present paper is sum- 
marized as follows: given model parameters, there is a 
maximum value of the gravitational potential produced 
by a star, above which no asymptotically de Sitter stellar 
configurations can be constructed. We show this both an- 
alytically and numerically. For example, the model with 
n = 1 and A « 2 gives $max ~ 0.1. This raises a warning 
sign for a class of f{R) theories, because neutron stars 
cannot be present in such gravity models. 

The underlying mechanism that hinders strong grav- 
itational fields around matter is explained essentially 
as follows :26]. Consider a static matter distribution. 
The Newtonian potential obeys the Poisson equation 
V^<I> r^ Gp, while the equation of motion for the scalar 
field implies V^x ^ Gp. From this, one can evaluate 
the excitation of the scalar degree of freedom around the 
matter distribution as 5x ~ ©($). If the de Sitter mini- 
mum is located very close to the point x = 1 , which cor- 
responds to i? = cx) in the effective potential, a slightly 
strong gravitational field will cause the problem of ap- 
pearance of a curvature singularity. 

Bearing the above evaluation in mind, let us comment 
on the other specific models of f{R) gravity. The model 
of Hu and Sawicki [20] and Starobinsky's one share the 
same structure of f{R) in the high-curvature regime, i.e., 
f{R) « i? - 2Aoff -I- C/R" with a > 0. Therefore, we ex- 
pect that the same problem arises in the Hu and Sawicki's 
model. The model by Appleby and Battye is character- 
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ized by [2l| 



fiR) 



R 

2" 



2a 



In [cosh(ai?) - tanh(6) sinh(ai?)] , (49) 



where a and b are parameters. Since x = df/dR = 
[1 + tanh(a_R— ^)]/2, a positive curvature singularity cor- 
responds to X = 1. Taking, for example, b = 1.5, we find 
that xi ~ 0.93 at the de Sitter minimum, which is very 
close to the dangerous curvature singularity (the result is 
independent of a). Since also in this model the effective 
potential is finite at i? = +00, we anticipate the same 
singularity problem. 

Our choice of the parameters in the present paper gave 
*i?inax ~ 0.1, for which neutron stars are unlikely to ex- 
ist. However, there still remains a possibility that more 
realistic stellar environments and matter profiles weaken 
the bound on the potential by a factor of 2 or 3, and at 
the same time make the chameleon mechanism work [36| . 
It is technically much more difficult to construct stellar 



configurations with a realistic equation of state, realistic 
energy densities, and realistic stellar environments. Such 
an elaborated modeling of relativistic stars might allow 
for <I>inax as large as, say, 0.3, but the parameter space of 
the theory will be very restricted. To conclude, f{R) the- 
ories that reproduce the correct behavior of weak gravity 
in the solar vicinity do not admit neutron star solutions 
without special care. 
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